Magnetic anisotropy of Co 2+ as signature of intrinsic ferromagnetism in ZnO:Co 
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We report on the magnetic properties of thoroughly-characterized Zni-xCo^O epitaxial thin films, 
with low Co concentration, x = 0.003 — 0.005. Magnetic and EPR measurements, combined with 
crystal field theory, reveal that isolated Co 2+ ions in ZnO possess a strong single ion anisotropy 
which leads to an "easy plane" ferromagnetic state when the ferromagnetic Co-Co interaction is 
considered. We suggest that the peculiarities of the magnetization process of this state can be 
viewed as a signature of intrinsic ferromagnetism in ZnO:Co materials. 

PACS numbers: 71.20.Be, 75.30.Gw, 76.30.Fc 



Spintronics, an emerging branch of micro- and na- 
noelectronics which manipulates the electron spin rather 
than its charge, has need for spin polarization compo- 
nents. In most spintronic devices, metallic ferromagne- 
tic (FM) materials are used to this end. However the 
physics of metal-semiconductor injection is incompatible 
with the concept of semiconductor devices, preventing 
their application [lj]. A suitable solution would be a FM 
semiconductor at room-temperature. 

The magnetic properties of diluted magnetic semi- 
conductors are due to the substitution of cations by 
transition-metal (TM) ions, and have been extensively 
studied for at least five decades Co-doped ZnO -a 
possible candidate for high-T c FM semiconductors - has 
attracted much interest from both theoretical and expe- 
rimental points of view. Yet, there is an ongoing debate 
about its magnetic properties. Early theoretical studies 
using the local spin density approximation (LSDA) for 
Zni-zCoaO found it to be a FM semimetal [3J. Contrary 
to this, more recent LSDA calculations [UTa] on large 
supercells detected a competition between FM and anti- 
ferromagnetic (AFM) interactions, i.e. an AFM or spin- 
glass groundstate. 

Experimentally, high-T c FM phases in Zn 1 _ 2 .Co :c O 
(x = 0.1-0.25) were found in thin films produced by 
pulsed laser deposition Q, by the sol-gel method 0], 
and by rf magnetron co-sputtering Q. They were also 
found in bulk single crystals prepared by implantation 
9]. Controversially, AFM correlations between TM ions 
and the absence of any FM bulk phases were observed 
in Zni-^Co^O (x — 0.005-0.15,0.2) samples fabricated 
by precursor decomposition in polycrystalline pow- 
der samples as well as in thin films fill ]. 

In this rather contradictory situation a major question 
which arises is whether a reliable identification of an in- 
trinsic FM phase of ZnO doped by Co is possible at all. 

Here we address this question on both experimental 



and theoretical grounds. We argue that such an identifi- 
cation requires a thorough examination of the magnetic 
properties of Co 2+ ions in the ZnO lattice, and in par- 
ticular the magnetic anisotropy of cobalt. By EPR and 
magnetic measurements, we first prove that Co 2+ , which 
has a spin S = 3/2, shows a huge single ion anisotropy of 
DS Z type, with D = 2.76 cm -1 . We then validate this 
result theoretically by combining crystal field theory with 
an estimate of the crystal field parameters. Theory and 
experiment clearly demonstrate that Co substitutes Zn 
in our samples. Finally, using a simple model, we show 
that a FM ZnO:Co would be an "easy plane" ferroma- 
gnet exhibiting a peculiar magnetization process which 
offers a simple and reliable way to identify the intrinsic 
FM phase of ZnO:Co. 

We focus on the magnetic anisotropy of isolated cobalt 
in ZnO and present details of the magnetic properties 
of epitaxial thin films with very low Co concentration 
varying from x=0.003 to x=0.005 \v^. The lfim thick 
samples were grown on sapphire substrate by plasma- 
assisted MBE. 2D growth is achieved for growth tempe- 
rature of 560°C, i.e. 50°C higher than the optimal growth 
temperature used for ZnO, resulting in streaky RHEED 
patterns. For this range of Co composition, the rocking 
curve FWHMs are in the range of oj ~ 0.15° along (002), 
(-105) and (105). The low ui values measured both for 
(-105) and (105), as well as their similarity, indicate a 
large column diameter, close to l^tm. The c-axis of the 
wurtzite structure is perpendicular to the film plane. The 
conductivity of the films is n-type, with residual carrier 
concentrations n e < 10 18 cm~ 3 , a doping level well below 
the Mott transition. 

Neglecting the hyperfine interaction, the 4 A2 ground- 
state of Co 2+ at a tetrahedral site of the ZnO host lattice 
is described by the following S — 3/2 spin Hamiltonian 



H S pin = ^Bg\\H z S z + tiB9±{H x S x + HySy) + DS 2 . (1) 
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The magnetic state of Co 2+ can thus be parameterized 
by only three constants : the two ^-factors, g\\ (H \\ c) 
and g± (H _L c), and the zero-field splitting constant 
D. Henceforth, we will use (Q]) with the values inferred 
from our experiments (see below), namely gu — 2.236, 
g± = 2.277 and D = 2.76 cm -1 , to compute (i) the ma- 
gnetization and the magnetic susceptibility of ZnO:Co in 
a simple statistical model of an ensemble of independent 
spins ; (ii) the angular and field dependences of the EPR 
spectra; (iii) the magnetization of small FM clusters. 
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Fig. 1: a) Field dependence of magnetization for a 
Zn^Coi-zO sample with x=0.0028 at T=2K. b) Inverse of 
the magnetic susceptibility for the same sample vs tempera- 
ture at /ioH=l T ; the inset shows x(T). Full and open circles 
are experimental data for H _L c and H \\ c, respectively. Solid 
lines are computed according to the model discussed in the 
text. 

First, we discuss the results of the magnetic measure- 
ments of thin ZnO:Co films. The substrate was cut into 
thin rectangular platelets of 3 by 3 mm 2 . In each expe- 
riment, up to 10 platelets were piled up to increase the 
paramagnetic component of the signal. In addition, un- 
doped ZnO films, deposited on the same sapphire sub- 
strate, were examined. The latter films served as a re- 
ference of the diamagnetic contribution of the measured 
signal. Measurements were performed using a Quantum 
Design MPMS XL magnetometer between 300K and 2K 
in magnetic fields up to 5 T. The field dependence of the 
magnetization taken at 2K and the temperature depen- 



dence of the inverse of the magnetic susceptibility measu- 
red at 1 T for a Zn^Coi-zO sample with x = 0.0028, for 
two orientations of the magnetic field, H _L c and H || c, 
are shown in Figs, la and lb, respectively. 

As is clear from Fig. 1, both the magnetization and the 
susceptibility curves reveal a significant magnetic aniso- 
tropy. This anisotropy can be referred to as being of an 
"easy plane" type ; i.e. for a given magnetic field H, the 
magnetization M _L c is greater than M || c. As expected, 
the susceptibility curves show a paramagnetic behavior. 
Their deviation from a Curie law finds a natural expla- 
nation in a S = 3/2 model with a single ion anisotropy, 
where the only adjustable parameter is the concentration 
of Co 2+ ions. 
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Fig. 2: Angular dependence of the effective g-factor of Co 2+ . 
9 = corresponds to H || c. The open circles are the experi- 
mental data ; the solid line is a fit based on Eq[Tj 



We now turn to the results of low-frequency EPR. An 
EMX Bruker spectrometer was used to collect spectra 
in the X-band (y = 9.4 GHz) and in the temperature 
range 4—300 K. A single line with a partly resolved hf 
structure for H || c was observed below ~ 100 K. This 
clearly indicates that D ^> hv. As the temperature is lo- 
wered to 4 K, the line intensity increases monotonically, 
and roughly follows a simple Curie law, ~ 1/T, indica- 
ting that the observed EPR signal is due to the low-lying 
doublet S 2 = ±1/2 of a S= 3/2 ground-state manifold 
14j. Note that the latter observation also allows us to 
determine the sign of D, D > 0. 

In Fig. 2 the angular dependence of the apparent (ef- 
fective) g-factor of Co 2+ in the ZnO lattice is shown, 
being the angle between the c-axis and the applied field 
H. The extracted values, g» — 2.236 and g e s± = 4.554, 
are very close to those obtained previously for single crys- 
tals [3] and thin films [l5| of ZnO:Co. The measured 
5off_L can be assigned to the lowest Co 2+ doublet, yielding 
g cB± ~ 2g x (l - (3/64) (^/L>) 2 ) for D > hv [lf| , which 
with reasonable accuracy reduces to g± = g c s±/2- The 
above results also suggest that the Zeeman part of the 
Co 2+ spin Hamiltonian is practically isotropic (pii ~ g±), 
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and hence cannot be responsible for the magnetic aniso- 
tropy of ZnO:Co. We therefore need more information 
about the zero-field splitting constant D which, accor- 
ding to the temperature dependence of the X-band EPR 
signal, can be estimated to 2D — 5.5 ± 0.3 cm" 
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Fig. 3: a) Transmission EPR spectra of ZnO:Co taken at 
A = 305/xm and 8 = 63°. b) Energy-level splitting of Co 2+ in 
magnetic field at 9 = 63° computed according to Eq[T] The 
vertical arrows indicate the observed EPR transitions. 



We have undertaken high-frequency EPR measure- 
ments on Zni-^Co^O (x = 0.01 — 0.02) single crystals 
in a large range of wavelengths A = 3mm — 0.3mm. 
Magnetic fields up to 35 T were provided by the pul- 
sed field facility in Toulouse. Representative transmission 
EPR spectra of ZnO:Co taken at the shortest wavelength 
(A = 305^im) are displayed in Fig. 3, together with the 
computed energy level splitting. The observation of the 
"forbidden" transitions with Am = ±2 and Am = ±3 
confirms the presence of a large zero-field splitting term 
in the spin Hamiltonian. Most importantly, these experi- 
ments allowed us to measure directly the zero-field split- 
ting constant, D = 2.76 ± 0.01cm -1 . 

To understand the origin of the magnetic anisotropy 
of Co 2+ , we use and extend the standard crystal field 
(CF) theory 13j, [l7| and estimate the CF parameters. 
The Co 2+ ion is in the d 7 configuration and Hund's rule 
coupling is sufficiently large to restrict the states essen- 
tially to S — 3/2, i.e. the A F and 4 P states. These states 
split due to the CF potential i?cF and interact by spin- 
orbit coupling Hso = XL ■ S : 



H = H 



Coul 



H, 



CF 



■Hi 



so 



Hqf — H cu b + H 



trig 



(2) 



The 4 P states are 15B higher in energy (B being the 
Racah parameter) . The trigonal symmetry of the CF re- 
quires three parameters which, as in Refs. dzj] and [H, 
we denote A, v, and v'. 

The importance of the non-diagonal matrix element 
v' was first pointed out by MacFarlane [l~7| but was not 



thoroughly treated in standard text books 13[ . MacFar- 
lane's formula 17f is a good approximation for strong 
cubic crystal fields (A ^> 15B, as for Cr 3+ complexes) 
but is not applicable in the present case where 15B 3> 
A ^> v, v', A. Here, in order to calculate <?|| , g± and D, we 
derive a perturbative formula for the parameters of the 
effective Hamiltonian (g s = 2.002, k : reduction factor) 
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Eq. ([3]) is more concise and hence more practicable than 
its counterpart derived in 

Numerical diagonalization of the Hamiltonian (|2J) wi- 
thin the subspace of the 4 F and 4 P states (a 40 x 40 
matrix) yields the results shown in Table fl] For the CF 
parameters obtained from optical measurements [l8[ , the 
anisotropy and gyromagnetic factors agree very well with 
our measured values. For the same parameters the per- 
turbative results of Eq. ^ are = 2.27, g±^ = 2.30, and 
2D = 4.34 cm -1 , and thus in good agreement with the 
numerical diagonalization. 
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Tab. I: Measured EPR data, compared to those calculated 
from CF theory with parameters from optics [Isl ] or esti- 
mated from the Harrison approach. In all calculations the 
Racah parameter B = 760 cm -1 and the spin-orbit coupling 
A = —143.3 cm -1 [ll| were used. The energy unit is inverse 
centimeters. 

To obtain more microscopic information about the CF 
parameters A, v and v' , we consider the hybridization 
contribution to the d- level splitting E m = (tprf m ) 2 /A p( i, 
with the Slater-Koster hopping parameter t„dm (fn = a 
or 7r) parameterized according to Harrison [20j, [2lJ . Ad- 
justing the charge-transfer energy A p d to the optical cu- 
bic splitting A gives A p d = 3.6 eV and k = 0.85 and 
leads to reasonable values of g\\ , g±, and 2D (last column 
of Table U). So, we may give to the CF parameters that 
were derived from the optical measurements a microsco- 
pic foundation. Also, these parameters are very specific 
to the tetrahedral environment of Co with trigonal split- 
ting which suggests that our experimental data can only 
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be explained if Co replaces Zn 
3.5 
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Fig. 4: Magnetization of FM Co 2+ clusters (N = 2, 4, 6) as a 
function of the applied magnetic field, calculated at T = 0.1K 
for J = -21cm -1 . 

We finally discuss the possible magnetic states that 
would arise if the individual Co 2+ ions in the ZnO lat- 
tice coupled ferromagnetically. This can be studied by 
adding a Heisenberg term JSiSj with J < 0) to the 
Hamiltonian ((T|) [22J. The magnetization curves, obtai- 
ned by exact numerical diagonalization of small S = 3/2 
clusters (with N — 2,4,6), shown in Fig. 4, are rather 
insensitive to the exact value of the FM coupling, provi- 
ded that — J ^> D. As can be seen from the figure, M(H) 
strongly depends on the orientation of the magnetic field : 
for H _L c, saturation is reached at very low fields; for 
H || c, by contrast, the magnetization rises at first essen- 
tially linearly (for N 3> 1) and saturates only at a critical 
field, H c . In the limit of N >■ 1, we have g\\iiBH c = 2D, 
and hence ^i Q H c = 5.3 T for FM ZnO:Co. We also note 
that the saturated magnetization does not depend on the 
magnetic field orientation. To our knowledge the above 
described magnetic behavior has not been seen yet. In 
contrast to this prediction an "easy axis" magnetic ani- 
sotropy 



23, 24 



251 ] or the absence of any anisotropy [6j 
were found in ferromagnetic polycrystalline ZnO :Co thin 
films, clearly indicating that Co electronic states other 
that 3d 7 are responsible for the observed ferromagnetic 
signal. 

In summary, we have presented an exhaustive study 
of the magnetic properties of Co-doped ZnO thin films 
in a low-concentration regime. Our experimental results 
and theoretical calculations clearly demonstrate a strong 
anisotropy of Co 2+ ions in the ZnO lattice which affects 
the magnetic ground-state of ZnO: Co leading to an "easy 
plane" ferromagnet. In the absence of consensus regar- 
ding the magnetic properties of ZnO: Co we argue that 
the study of its magnetic anisotropy offers simple criteria, 



both experimental and theoretical, for the identification 
of the intrinsic ferromagnetism in this material. 
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